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The application of Type la supernovae (SNe la) as distance indicators in cosmology calls for 
a sound understanding of these objects. Recent years have seen a brisk development of astro- 
physical models which explain SNe la as thermonuclear explosions of white dwarf stars. While 
the evolution of the progenitor is still uncertain, the explosion mechanism certainly involves the 
propagation of a thermonuclear flame through the white dwarf star Three-dimensional hydro- 
dynamical simulations allowed to study a wide variety of possibiUties involving subsonic flame 
propagation (deflagrations), flames accelerated by turbulence, and supersonic detonations. These 
possibilities lead to a variety of scenarios. I review the currently discussed approaches and present 
some recent results from simulations of the turbulent deflagration model and the delayed detona- 
tion model. 
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1. Introduction 

Thermonuclear supernova explosions are an astrophysical model for the astronomical class of 
Type la supernovae (SNe la henceforth). These objects are of interest in many fields of astrophysics 
and astronomy. Being one of the main sources of iron group elements, SNe la contribute to the 
chemical evolution of galaxies (e.g., [1]). They affect star formation and drive shock waves in the 
interstellar and intergalactic media. 

The main driver of SN la research over the past years, however, has been their application 
in observational cosmology. Here, SNe la were employed as distance indicators (as put forward 
by [2]). At redshifts above 0.5, a significant deviation from the linear Hubble law was noticed 
which led to the spectacular interpretation of the Universe currently undergoing an accelerated 
expansion [3, 4]. The determination of the force driving this acceleration is perhaps one of the 
greatest challenges in contemporary physics. For the time being, it is parametrized as "dark energy" 
(e.g., [5]). The simplest form of dark energy is a cosmological constant, but more complicated 
contributions to the energy-momentum tensor in the Einstein equations are also conceivable. A 
first step to determine the nature of dark energy would be to constrain its equation of state. SNe la 
seem to be a suitable tool for this task and currently two major campaigns [6, 7] apply them in 
distance determinations of hundreds of supemovae out to redshifts of z ~ 1 . The large number of 
observations is necessary to reduce the statistical errors because putting tight constraints on dark 
energy equation of state requires a high accuracy of the distance determinations. 

From a theoretical point of view, however, the applicability of SNe la as distance indicators is 
still not satisfactorily answered. SNe la are remarkably uniform events by astrophysical standards, 
but evidently no standard candles. Only a calibration of the distance measurements according to 
empirical correlations between observables provides the necessary accuracy for the determination 
of cosmological parameters. Such calibrations may be afflicted with systematic errors. Being 
derived from a set of nearby well-observed SNe la, there is no guarantee that they perform well 
for supernovae at high redshifts, too. Getting a handle on these uncertainties is one of the goals of 
modeling SNe la. 

2. Astrophysical modeling 

The cornerstones of the astrophysical model of SNe la derive from the fundamental character- 
istics of these events: 

• Evidently, SNe la belong to the most energetic cosmic explosions, releasing about 10^^ erg 
of energy. For a short period of time they can outshine an entire galaxy consisting of tens of 
billions of stars. 

• SNe la spectra are characterized by lacking indications of hydrogen and helium which to- 
gether with a pronounced P Cygni silicon line at maximum light classifies these objects [8]. 
Lines of intermediate-mass elements (IMF, such as Si, Ca, Mg and S) and oxygen are ob- 
served in near-maximum light spectra (e.g., [9]). 

• SNe la form a class of remarkable homogeneity with respect to observed lightcurves and 
spectra (e.g., [2]). 
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Assuming supernovae to generally originate from single stellar objects, only their gravitational 
binding energy, released in a collapse towards a compact object [10], or its nuclear energy, released 
in explosive reactions [11], come into consideration as possible energy sources. In the particular 
case of SNe la, no compact object is found in the remnants excluding the first possibility. The 
homogeneity of the class of SNe la and the fact that no hydrogen is found in their spectra provides 
a strong hint that the object undergoing the nuclear explosion may be a white dwarf (WD) star 
consisting of carbon and oxygen (C+0). 

Lightcurves of SNe la rise over a time scale of several days and decline over months. It is 
therefore clear that they cannot be powered directly by the explosion since the temperatures fall off 
much too rapidly in the expansion. This problem was solved by noting that the ^^Ni produced in 
large amounts in the explosive thermonuclear burning provides the energy source for the optical 
event by radioactive decay to ^^Co and ^^Fe [12, 13]. 

2.1 Progenitor evolution and ignition 

A single WD is an inert object. How can it reach an explosive state? The only way to introduce 
the necessary dynamics into the system is to assume it to be part of a binary system and to gain 
matter from the companion. Several models have been proposed for this progenitor evolution. 
Here, we will discuss only models resulting from the so-called single-degenerate Chandrasekhar- 
mass scenario, which has received by far most attention recently. For a summary of alternative 
models, we refer to [14]. 

In the single degenerate scenario [15, 16, 17], the WD accretes matter from a non-degenerate 
companion (either a main sequence star, or an AGB star). This idea was recently supported by the 
detection of the potential companion of Tycho Brahe's supernova of 1572, which is a solar-type 
star [18]. By accreting material from this companion and burning it hydrostatically at the surface, 
the WD may reach the Chandrasekhar mass Mch- Limiting the fuel available in the explosion to 
Mch ~ 1.4Mo, this Chandrasekhar-mass model appears particularly favorable since it provides a 
natural explanation for the striking uniformity of SNe la in the gross observational features. At the 
same time, it is afflicted with great uncertainties. Achieving a stable mass transfer in the progenitor 
binary system to build up a Chandrasekhar mass WD is highly non-trivial (e.g., [19]) and the 
observational evidence for such systems is sparse. 

When the WD approaches the Chandrasekhar limit, the density at the center of the WD in- 
creases rapidly so that fusion of carbon ignites. Contrary to the situation in main sequence stars, 
the degenerate material of the WD does not allow for moderation of the burning by expansion. 
Heat transport is achieved here by convection giving rise to a stage of convective carbon burning 
that lasts for several hundred years. This phase is terminated by one or more small spatial regions 
undergoing a thermonuclear runaway, marking the birth of a thermonuclear flame and the onset 
of the explosion. The convective burning stage and the conditions at flame ignition are extremely 
hard to model both analytically and numerically. Therefore the exact shape and location of the first 
flame spark(s) is not yet well constrained [20, 21, 22, 23, 24]. Does ignition occur in a single spot 
or in multiple sparks with a stochastic distribution? Does it take place at the center of the WD or 
do pre-ignition convective motions lead to large asymmetries? Evidently, the ignition structure is a 
crucial initial parameter in multi-dimensional explosion models [25, 26, 27, 28, 29, 30]. 
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2.2 Flame propagation and explosion 

The goal of SN la explosion models is to follow the propagation of the thermonuclear flame 
from its ignition near the center of the WD outwards and to determine the nuclear energy release 
and the structure of the ejected material. 

Hydrodynamics allows for two distinct modes of flame propagation. One is the subsonic 
deflagration in which the flame is mediated by the thermal conduction of the degenerate electron 
gas and the other is a supersonic detonation in which the burning front is driven by shock waves. 
Either one of these modes or a combination of both have been suggested in different explosion 
models: 

• The prompt detonation model [31] attempts to explain SNe la by a detonation ignited at the 
center of the WD and propagating outward. This produces enough energy for a SN la event. 
However, ahead of a supersonic detonation wave, the fuel cannot expand and is therefore 
incinerated at the high densities of an equilibrium white dwarf. This results in the almost 
complete conversion of the material to nickel-peaked nuclear statistical equilibrium [31], 
which is in conflict with the intermediate mass elements observed in SN la spectra. Such 
nucleosynthetic problems rule out pure detonations as a standard model for SN la explosions. 

• The deflagration model [32] assumes the flame to propagate in the subsonic deflagration 
mode. The laminar burning speed of the deflagration flame is determined by microphysical 
transport processes. For conditions of carbon burning in C+0 WDs it is highly subsonic [33] 
and therefore the flame propagates far too slowly to explain SN la explosions. The expansion 
of the star quenches burning before the WD gets unbound. On the other hand, this model can 
cure the problem of nucleosynthesis, since rarefaction waves travel ahead of the flame with 
sound speed and lower the fuel density prior to burning. Thus, the material can partly be 
processed into intermediate mass elements. The deflagration model undergoes a significant 
improvement when multidimensional eff'ects are taken into account as will be discussed in 
detail in Sect. |. 

• The delayed detonation model [34, 35, 36] unites the advantages of the deflagration and the 
detonation models. Burning starts out in the slow deflagration mode pre-expanding the star. 
At some point, a transition from the initial subsonic deflagration to a supersonic detonation 
takes place. This detonation is an easy way to explain the energy release necessary for a 
SN la explosion. The important notion in this model is that a detonation in low density 
fuel (pre-expanded in the deflagration stage) can lead to only partial burning and is therefore 
capable of generating intermediate mass elements. A detailed account of this model will be 
given in Sect. |5[ 

• The pulsational delayed detonation model [37] is similar to the delayed detonation model in 
the sense that it combines an initial deflagration with a later detonation. The flame is assumed 
to propagate in the initial deflagration phase with its laminar burning speed and pre-expands 
the star. Due to the slow flame velocity, the burning front stalls and fails to unbind the star. 
The WD then re-contracts giving the interface between burnt and unburnt material enough 
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time to mix and to become nearly isothermal. Compressional heating finally triggers a deto- 
nation at densities that are lower than that prior to the first expansion phase. The assumption 
of the flame propagating with the pure laminar burning velocity in the deflagration phase 
originally made in this model seems unrealistic because of the flame instabilities and the re- 
sulting turbulent flame acceleration. Recent multidimensional deflagration models [38, 39] 
demonstrated that taking these effects into account, the star is likely to get unbound instead 
of recontracting. 

3. Numerical modeling 

Numerical models of SN la explosions have to face three major challenges. The vast range 
of relevant length scales cannot be resolved in computational models in the foreseeable future. It 
extends from the radius of the WD star (~2000km at the onset of the explosion and expanding in 
the process) down to the flame width which is well below one centimeter and, including turbulence 
effects, to the Kolmogorov scale of less than a millimeter. Apart from this problem, numerical 
models need to take into account inherently three-dimensional physical phenomena and to solve 
the equations of nuclear burning. To meet all these requirements in a single simulation will be im- 
possible in the foreseeable future. Therefore the problem has been tackled in different approaches. 

The first path towards SN la explosion modeling is to restrict the simulations to only one 
spatial dimension. Here, in principle a resolution of the relevant scales is achievable and a de- 
tailed description of the nuclear reactions is feasible. However, this approach fails to consistently 
incorporate crucial three-dimensional physical mechanisms. 

In multi-dimensional simulations (see [40, 41] for recent reviews), the computational costs of 
modeling the explosion hydrodynamics is prohibitive to directly resolve all relevant scales as well 
as details of the nuclear processes. While the latter may be improved in the forthcoming years, 
and is meanwhile separated from the actual explosion simulations still maintaining a reasonable 
accuracy (for an approach based on reconstructing the nucleosynthesis from tracer particles, see 
[42, 43]), even a drastic increase in computational capabilities will not allow for a resolution of 
all relevant scales in multi-dimensional simulations. Consequently, additional modeling effort is 
required in order to implement a consistent description of flame propagation in these simulations 
(see Sect. |). 

A third approach is to study specific effects on a limited range of spatial scales, in order to 
validate assumptions and improve modeling techniques of the large-scale SN la simulations. 

3.1 Objectives of modeling Type la supernovae 

Three-dimensional numerical models of Type la supemovae strive for a self-consistent astro- 
physical description of these astronomical events. The goal is to achieve a viable model by starting 
out from fundamental physical laws ("first principles") and to formulate the model without intro- 
ducing free parameters. Avoiding tunable parameters allows the simulations to gain a high level 
of predictive power. Only this way, a direct comparison of synthetic observables derived from 
the simulations with actual SN la observations facilitates a thorough validation of the modeling 
approach. 
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This aim of multidimensional models is significandy different from earlier one-dimensional 
modeling approaches. Here, the flame propagation velocity, which is determined by inherently 
multidimensional eff'ects such as instabilities and turbulences, cannot be consistently determined. 
It thus introduces a free parameter which has significant impact on the result of the simulated explo- 
sion. With such a powerful tunable parameter an impressive agreement between one-dimensional 
simulations and observations could be achieved. This in turn tells us a great deal about the global 
properties of the explosion mechanism. Successful one-dimensional models (e.g. W7 [44]) bring 
forward the average flame velocity that needs to be attained in different stages of the explosion and 
remain the benchmark of the (spherically averaged) chemical composition of the explosion ejecta. 

Successful multidimensional models, however, would diverge from this way of tackling SN la 
modeling as an inverse problem. Modeling the explosion from first principles, they make predic- 
tions of observables. This paves the way to investigating questions related to the cosmological 
use of SNe la: What is the precision of SNe la as distance indicators? What is the reason for 
the SN la diversity? How does the diversity relate to variations in the properties of the progeni- 
tor system? How do these variations translate (via the explosion process) to observable features? 
How can one improve the calibration techniques of SN la cosmology? With ongoing observa- 
tional campaigns, projects currently starting and future surveys, the number of observed SNe la 
will increase from currently a few hundred per year to several thousands. This calls for a sound 
theoretical understanding of SNe la as an astrophysical phenomenon. A rapid development of the 
field of multidimensional SN la simulations gives hope that these models will be able to tackle the 
questions raised above in the not too far future. 

3.2 Current status of SN la modeling 

The zoo of variants of the Chandrasekhar-mass single-degenerate SN la model recently studied 
in multidimensional simulations comprizes 

• the turbulent deflagration model (see also Sect. ^), which is so-far best studied in multidi- 
mensional simulations [38, 45, 39, 46, 47, 25, 26, 28, 48]. The advantage of this model is that 
the explosion process can be formulated free of tunable parameters. Therefore it is possible 
to compare the outcome of simulations directly to observations in order to assess the validity 
of the model. Systematic tests of the initial parameters of the model on the outcome of the 
exploson process [49, 43] and their implication for SN la cosmology [50] are possible. 

• the delayed detonation model [51, 52, 53, 54], which mainly suffers from the uncertain 
mechanism providing a deflagration-to-detonation transition in thermonuclear supernovae. 
Fixing this uncertain parameter to a physically motivated hypothesis, however, the results of 
corresponding numerical simulations look promising (see Sect. ||). 

• the gravitationally conflned detonation model, which arises from asymmetric deflagration 
flame ignitions. It has been claimed that in some cases the deflagration may fail to unbind 
the WD [27]. As a consequence of the off-center ignition, the ashes of the deflagration 
burning quickly float towards the surface of the WD. If still gravitationally bound, they start 
to sweep around the unburnt core of the star and collide on the far side. This collision has 
been claimed to trigger a detonation wave [27, 55, 56] (see, however, [28, 57]) propagating 
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inwards and burning out the core of the WD. Very energetic and bright events are expected 
from such a model. 

• the pulsational reverse detonation model, which follows the pulsational phase of the gravi- 
tationally bound WD if no detonation is triggered by the collision of ashes as in the gravita- 
tionally confined detonation scenario. Pulsational contraction has been suggested to trigger 
a detonation reviving the explosion in this case [58]. 

The goal of a fully self-consistent SN la model that agrees with all observatonal features is not 
reached yet. Some of the modeling approaches, however, seem very promising. The following part 
of the paper focuses on a detailed description of the first two of the above listed models. 

4. The turbulent deflagration model 

The turbulent deflagration model extends the laminar deflagration model by considering ef- 
fects of instabilities and turbulence on the flame propagation. These lead to a significant accelera- 
tion of the burning and to viable explosions of the WD star. 

4.1 Turbulent deflagratlion in thermonuclear supernovae 

The major effect accelerating the flame is due to the buoyancy unstable flame propagation 
from the center of the star outwards. It leaves behind light and hot ashes below the dense fuel - a 
density statification inverse to the gravitational acceleration. In its non-linear stage, the Rayleigh- 
Taylor instability leads to the formation of mushroom-shaped burning bubbles raising into the 
fuel. The Reynolds number typical for this situation is as high as 10^^. Clearly, shear (Kelvin- 
Helmholtz) instabilities at the interfaces of these bubbles will generate turbulent eddies which then 
decay to smaller scales forming a turbulent energy cascade. The flame will interact with these 
eddies down to the Gibson-scale at which the turbulent velocity fluctuations become comparable 
to the laminar flame speed. Below the Gibson scale, the flame burns faster through turbulent 
eddies than they can deform it, and the flame propagation is thus unaffected by turbulence there. 
This interaction corrugates the flame increasing its surface, enhancing the net burning rate, and, 
consequently, accelerating the effective flame propagation speed. 

In the numerical implementation of the deflagration model, the relevance of turbulent effects 
amplifies the scale problem since the turbulent cascade extends to the extremely small Kolmogorov 
scale (< 1 mm) where the turbulent energy is dissipated into heat. The flame interaction with the 
turbulent cascade down to the Gibson scale must be taken into account. Current 3D simulations 
capturing the entire star reach resolutions around one kilometer while the Gibson scale is of the 
order of lO'^cm at the beginning of the explosion and decreases steadily. For large-scale multi- 
dimensional SN la simulations this has three consequences: 

1. The internal flame structure cannot be resolved. Thus, an effective flame model has to be 
applied and complementary small scale simulations are required. 

2. It is not possible to fully resolve the interaction of the flame with turbulence. Therefore 
modeling of the effects on unresolved scales is necessary. 
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3. Assumptions about the flame properties on unresolved scales [59, 49, 60, 61, 62] have to be 
validated in separate small-scale simulations. 

Multi-dimensional models are inevitable to consistently determine the turbulent flame propagation 
velocity. Given the wide range of scales on which the flame is affected by turbulence, this is an 
ambitious project, additionally challenged by the lack of resolution of the thermonuclear flame 
structure. For both problems, different approaches have been taken, guided by the theory of turbu- 
lent combustion in terrestrial flames (see [63]). 

4.2 Flame model 

Two major strategies to tackle the problem of the unresolved internal flame structure can be 
distinguished. On the one hand, a flame capturing technique [64] mimics flame propagation by an 
artificial diffusion mechanism which broadens the internal flame structure to a certain number of 
computational grid cells. It has been applied in SN la explosion simulations [39, 51, 65, 27, 56]. 

On the other hand, a completely different approach [66] treats the flame as a sharp discontinu- 
ity separating the fuel from the ashes. It is numerically represented applying the level-set technique 
[67]. Here, the flame front is associated with the zero level set of a scalar function G represent- 
ing the distance from the interface. A model for flame propagation based on this technique [68] 
was modified for thermonuclear flames in SN la explosions [66, 69]. This scheme was applied in a 
number of simulations [45, 38, 70, 46, 47, 71, 72, 26, 43, 25, 28, 48], and the simulations presented 
below are based on it. 

4.3 Turbulent combustion model 

The wide range of scales involved in turbulent combustion in thermonuclear supernovae makes 
direct simulations virtually impossible. Only parts of the flame/turbulence interaction range and of 
the resulting flame surface enlargement can be resolved on the computational grid. This deficit is 
usually compensated by attributing an effective turbulent flame speed Si to the unresolved flame 
front, which must be determined by theoretical considerations. 

One of the cornerstones of the theoretical description of turbulent combustion is the notion 
of different regimes of flame/turbulence interaction [73]. These regimes are distinguished by the 
ability of turbulent eddies to penetrate the internal flame structure. Since the Gibson scale is much 
larger than the flame width for most parts of the SN la explosion, this will not be the case here 
and accordingly the combustion falls into the regime of wrinkled and corrugated flamelets (e.g., 
[63]). In this regime, the full flame structure is corrugated by the interaction with turbulence and 
the resulting surface enlargement accelerates its propagation. The flame propagation completely 
decouples from the microphysics of the burning for sufficiently strong turbulence [74]. It is entirely 
determined by the turbulent velocity fluctuations, that is, the effective turbulent flame speed ^t is 
proportional to the turbulent velocity fluctuations v'. 

One of the challenges of deflagration models of SN la explosions is thus to determine these 
velocity fluctuations correctly. Since the resolution in multi-dimensional simulations is insufficient 
to resolve the phenomena directly, modeling approaches have to be taken. The assumption of the 
flame propagation being entirely driven by buoyancy instabilities [39, 65, 27, 56] falls short of 
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reproducing the effects of a turbulent cascade which will develop due to large-scale shear insta- 
bilities and dominate the turbulence properties at small scales. Such effects are taken into accunt 
only with appropriate subgrid-scale turbulence models [75, 72]. Guided by the technique of Large 
Eddy Simulations (LES), this model determines the turbulence energy on unresolved scales based 
on conservation laws. It is applied in the simulations discussed below. 

In the very late stages of the SN la explosion, the fuel density drops due to the expansion of 
the WD to values where the flame width becomes broader than the Gibson length. Then, turbulence 
penetrates the internal structure of the flame and it enters the regime of distributed burning [73]. 
Including this burning stage into thermonuclear supernova simulations [46, 76] affects the latest 
stages of deflagration burning. 

4.4 Example: a highly resolved simulation 

Several simulations based on the implementation described above, both in two and three spatial 
dimensions, have been presented [77, 45, 38]. In the 2D simulations, numerical convergence in the 
global quantities was demonstrated. For the implementation on a co-expanding computational grid, 
a similar result was found [71]. The numerical convergence naturally arises from the interplay of 
the resolved flame front representation with the turbulent subgrid-scale model. Ideally, a lack of 
resolution of large-scale features in the flame front representation should be compensated by an 
increased turbulent flame propagation velocity determined from the subgrid-scale approach. Of 
course, a certain threshold of resolution will need to be exceeded in order to reach this regime in 
the numerical implementation. A consistent and reliable sub-grid scale modeling was the goal of a 
recent highly resolved simulation [48]. This model shall be described here in order to illustrate the 
typical flame evolution in the deflagration scenarios of thermonuclear supernovae. A full-star [47] 
multi-spot ignition [25] model was set up on a computational grid comprizing 1024 x 1024 cells. In 
combination with the nested moving grid approach [25] this facilitated an extremely fine-structured 
initial flame geometry (see the upper left panel of Figure |l|). 

Starting from this initial flame configuration, the evolution of the flame front in the explosion 
process is illustrated by snapshots of the G = isosurface at ? = 0.6s and t = 3.0s in Figure |l] 
(upper right and lower left panels). The development of the flame shape from ignition to ? = 0.6 s 
is characterized by the formation of the well-known "mushroom-like" structures resulting from 
buoyancy. During the flame evolution, inner structures of smaller scales catch up with the outer 
"mushrooms" and the initially separated structures merge forming a connected configuration (see 
snapshot at ? = 0.6 s of Fig. [I]). The continued development of substructure and the merger of 
features create a deflagration structure with a complex pattern. Burning and flotation drive the 
flame towards the surface of the WD. The fuel density drops as the flame moves outwards due 
to radial stratification and the overall expansion of the WD caused by the energy deposit from 
nuclear burning. Once the fuel density falls below the threshold for the production of intermediate- 
mass elements, nuclear burning ceases. This occurs first at the leading features of the flame and 
subsequently in more central flame regions. Finally, no burning takes place anymore. At this point 
the outer ash features have reached the surface layers of the ejecta (lower left panel of Fig. [T]). 
The following seconds in the evolution are characterized by hydrodynamical relaxation towards 
homologous expansion. This stage is reached to good approximation at about 10 s after ignition 
[71]. 
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Figure 1: Snapshots from a full-star SN la simulation starting from a multi-spot ignition scenario. The loga- 
rithm of the density is volume rendered indicating the extend of the WD star and the isosurface corresponds 
to the thermonuclear flame. The last snapshot corresponds to the end of the simulation and is not on scale 
with the earlier snapshots (from [48]). 

This three-dimensional evolution leads to a remnant of the explosion with characteristic prop- 
erties. The density structure has patterns from unstable and turbulent flame propagation imprinted 
on it (see the lower right panel of Fig. |l|), and ash regions extend to the outermost layers of the 
expanding cloud of gas. 

4.5 Comparison with observations 

Due to recent progress in deriving observables from multi-dimensional deflagration simula- 
tions, a direct comparison with details of observations of nearby SNe la has come into reach. 
Since the simulations contain no other parameters than the initial conditions, the question arises of 
whether they meet observational constraints. Such constraints result from the global characteris- 
tics derived from observations, observed lightcurves, and spectra taken from nearby SNe la. We 
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Figure 2: Bolometric lightcurve derived from the highly resolved simulation (black curves; solid is the 
"UVOIR-bolometric" lightcurve and the complete bolometric lightcurve is dotted). The blue dotted curves 
correspond to observed bolometric lightcurves [79] (from [48]). 

will present a comparison of the highly resolved simulation described above with observational 
expectations in the following. 

The global characteristics derived from SN la observations state that a valid explosion model 
should release -lO^'erg of energy and produce about 0.4 . . . 0.7 M© of ^^Ni in the nuclear burn- 
ing [78, 79]. However, there exists a large diversity in the observations ranging from the class of 
sub-luminous SNe la (like SN 1991bg with probably ~0.1Mo of ^^Ni) to super-luminous events 
(e.g. SN 199 IT with a ^^Ni mass close to IM©). The simulation under consideration here led to 
an asymptotic kinetic energy of the ejecta of ~8.1 x 10^" erg and produced ~ O.6M0 of iron group 
elements. It thus falls into the range of observational expectations, albeit on the weaker side. 

Lightcurves — Lightcurves of SNe la are sensitive to the energy release, the ^^Ni production, as 
well as to the distribution of elements in the ejecta. In Figure ^ a synthetic bolometric lightcurve 
derived from the above described simulation is compared with observed lightcurves [53]. It was 
calculated using the Stella code [80, 81, 82]. The synthetic bolometric lightcurve falls into the 
range of observed lightcurves of SNe la (indicated in Figure ^ by the extreme examples of the 
super-luminous SN 1991T and the sub-luminous SN 1991bg). With respect to brightness and 
shape it compares reasonably well with SN 1994D, a standard normal SN la. 

Spectra — A much harder test for the models is posed by the comparison of synthetic and ob- 
served spectra since these depend on details in the composition of the ejected material. A powerful 
diagnostic tool to compare SN la models with spectral observations is provided by the abundance 
tomography of SN 2002bo [83]. It makes use of spectra taken from this supernova with an ex- 
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Figure 3: Spherically averaged composition resulting from the hydrodynamical explosion simulation (solid 
lines) compared to the findings of the abundance tomography of SN 2002bo (dotted lines). Iron group 
element abundances are shown in red, intermediate-mass elements in green, and unburned material in blue 
(from [48]). 

traordinary good time coverage. Fitting this sequence of data with synthetic spectra unveils the 
composition of the ejecta in velocity space slice by slice (see Figure ||), since the photosphere 
moves gradually inwards with the expansion of the remnant. 

This abundance tomography of the ejecta can be compared with results of 3D models, when 
averaged over the angles (Figure ^). Qualitatively, the mixed composition of the ejecta found 
here is reproduced by deflagration SN la models in a natural way since these predict a distribution 
of burnt material within the rising bubbles. A problem was, however, that older models predicted 
large unbumt mass fractions in the central parts of the ejecta in disagreement with the observational 
results [84, 83]. The high-resolved simulation cures this problem by clearly reproducing the iron- 
group dominance in the low- velocity ejecta [48]. A good agreement between the simulation and 
the abundance tomography of SN 2002bo is therefore recovered in the inner parts of the ejecta. 
Above radii of ~ 10, 000km s~^ however, the simulation predicts the chemical composition to be 
dominated by unbumt material. This is in contradiction to observations which find the material still 
to be dominated by nuclear ashes, mainly intermediate-mass elements. It should however be noted 
that these outer regions of the ejecta contain only little mass (~0.25Mo) due to the low densities 
encountered here. Thus, the observed chemical composition is reproduced by the model for the 
larger part of the ejected material. 

4.6 Summary of the deflagration model 

The pure deflagration model has proven very successful in many respects. It provides the up 
to now only fully self-consistent thermonuclear supernova model. Apart from the initial conditions 
(flame ignition, WD structure and composition), which are expected to vary in nature, it describes 
the actual explosion process without tunable parameters. Perhaps the most striking success is that 
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it indeed leads to an explosion of the WD releasing energies that fall into the range of observational 
expectations. The structure of the ejecta seems to be partially consistent with compositions of SN la 
ejecta derived from spectral observations. At least the inner part of the ejecta looks similar to what 
is expected for the dimmer examples of normal SNe la. 

Despite these successes, there are shortcomings of the pure deflagration model. Varying the 
initial conditions in a range that seems physically plausible, it could not reproduce the full range 
of observed SN la diversity. In particular, the normal to bright examples of the observational 
SN la sample seem to be out of reach for the deflagration model. For these to be reproduced, 
^^Ni masses of the order of 0.7 M© need to be synthesized in the explosion. This could not be 
achieved thus far in the deflagration model. The second concern is with regard to the composition 
of the ejecta. The large-scale Rayleigh-Taylor instabilities lead to clumpy inhomogeneities in the 
structure of the explosion ejecta and to a strong mixing of different groups of reaction products 
in the angular averaged profile. While the results may still be consistent with observations in the 
central parts of the ejecta, the chemical composition of the outer parts seems to be at odds with 
the observation, where an intermediate-mass element dominance of the ejecta is detected out to 
large radii. It seems that the burning provided by the deflagration flame is insufficient, in particular 
in late phases of the explosion. Such an incomplete burning would explain the shortcomings in 
reproducing the outermost parts of the ejecta and the fact that the deflagration model fails to explain 
the brighter SN la events. Although the kinetic energy of the explosion ejecta falls into the range 
of observational expectations, it reaches only the lower end of what is expected for normal SNe la. 
Two interpretations are possible here: 

1. The deflagration model accounts for a certain peculiar sub-class of SNe la [85], namely 
rather dim events which show indications for unburnt material in their ejecta. 

2. The deflagration model is incomplete but provides an important building block for an ex- 
tended explosion model. This would explain why it is rather successful in reproducing the 
dimmer SN la events. For these, the deflagration phase may dominate the explosion process. 

In the following section, we elaborate on the second of these possibilities. 

5. The delayed detonation model 

A promising extension of the pure deflagration model is the delayed detonation model [35]. 
By assuming a transition of the flame propagation mode from subsonic deflagration to supersonic 
detonation in an advanced stage of the explosion process, it provides a way to accelerate and en- 
hance burning in the late phases. According to the discussion above, this may be a way to achieve 
a better agreement with SN la observations. Not only would it lead to a stronger explosion and 
to an enrichment of the outer layers of the ejecta with nuclear ashes, but it would also weaken the 
clumpy inhomogeneities in the ejecta by supersonically burning down the funnels in between the 
uprising Rayleigh-Taylor plumes. 

5.1 Deflagration-to-detonation transitions in SNe la 

The delayed detonation model, however, is still afflicted with severe uncertainties. Certainly 
the greatest hindrance for its success is that the physical mechanism providing a deflagration-to- 
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Figure 4: Turbulent velocity fluctuations u' at patches of the flame which have entered the distributed 
burning regime. The blue opaque, transparent, or wire-mesh surfaces correspond to the flame front. Volumes 
of high turbulent velocity fluctuations are rendered in red/orange. In the lower panels, the green arrow 
indicates the location of the maximum value of v' found in the simulation. For better visibility, white areas 
correspond to ash regions and fuel regions are shown in black in a plane intersecting with the maximum 
y'-value in the lower right panel (from [86]). 



detonation transition (DDT henceforth) in thermonuclear supernovae remains unclear [87]. Such 
transitions are observed in terrestrial combustion, but there they occur in the vicinity of obstacles or 
at the walls of the combustion vessel. Such boundaries do not exist in the astrophysical situation. 
Therefore, it seems reasonable to assume that a DDT occurs (if happening at all) in connection 
with the only significant change in the structure of the deflagration flame that occurs when it enters 
the distributed burning regime [88]. Remarkably, this takes place at fuel densities of the order 
lO^gcm"^, coinciding with the "transition density" that led to best results in parameterized one- 
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Figure 5: Fits to the histogram of the turbulent velocity fluctuations v' (from [86]). 



dimensional delayed detonation simulations. 

DDTs occurring in the distributed burning regime seem to be possible [89], but only if a 
very high turbulence intensity is retained in these late explosion processes. Turbulent velocity 
fluctuations of the order of 1000km s"^ seem to be necessary to trigger a DDT. While earlier two- 
dimensional simulations of the deflagration phase in thermonuclear supernovae failed to reach these 
values, a recent analysis of three-dimensional simulations [86] found that such high turbulent in- 
tensities are likely to be realized in late burning stages. An example is given in Figure ^ illustrating 
the situation in the highly resolved deflagration simulation described in Sect. ^. Patches of the 
flame front are color-coded which have entered the distributed burning regime and still feature 
high turbulent velocities. A histogram of the velocity fluctuations corresponding to the situation of 
Figure ^ is shown in Figure |5[ It features an extended high-velocity tail. An acceptable fit to this 
histogram is obtained with a probability density function following an exponential of a geometric 
Ansatz: 



P{v')^QMaQ{v'T+a2\. 



(1) 



For values of the fitting parameters ao, ai, and ai see [86]. From such an analysis it is evident that 
the necessary high turbulent velocity fluctuations are likely to be achieved at non-negligible patches 
of the flame front and a flame-driven DDT thus cannot be excluded. However, the microscopic 
mechanism of this transition remains largely unclear. A high value of the mean turbulent velocity is 
necessary, but certainly not sufficient [90]. Further studies of the microscopic properties of burning 
in the distributed regime are needed in order to settle the question of a DDT in thermonuclear 
supernovae. 
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Figure 6: Initiation and propagation of the detonation wave (white isosurface) in a delayed detonation model 
[53]. The deflagration flame is shown as blue isosurface and the extent of the star is indicated by the central 
plane mapping the logarithm of the density (from [53]). 

5.2 Delayed detonations in full-star simulations 

Despite the uncertainties of DDTs, an approach to address the question of delayed detonations 
in thermonuclear supernovae is to assume such a transition and to artificially trigger the detonation 
wave in large simulations. The question to be addressed by this approach is whether the outcome 
of such models is consistent with observations. Prescribing the DDT to occur at an arbitrary instant 
and a location near the center led to promising results [51]. However, such an approach suffers 
from the arbitrariness of the assumed DDT and has little predictive power. 

A way to assess the potential of the delayed detonation model is to fix the DDT to some 
physically motivated guess [53]. According to the discussion above, one possibility is to assume 
the DDT to occur at the location where the flame first enters the distributed burning regime [52, 53]. 
Although this is a very simplifying assumption, it fixes the unknown DDT parameter and allows 
to compare the results of the simulations with observational expectations. This was done both in 
two-dimensional [52] and three-dimensional [53] simulations. Particular care was taken to prevent 
the detonation from crossing ashes left behind by the previous deflagration stage. This requirement 
follows from a recent analysis of the propagation of detonations in WD matter [91]. It can easily 
be implemented by following the detonation wave with a level-set approach [52, 28, 92]. 

A typical evolution of a three-dimensional full-star simulation is shown in Figure ^. After 
ignition, the deflagration phase proceeds in a way similar to the pure deflagration model. The 
burning bubbles grow by flame propagation and rise buoyantly towards the surface of the WD. Due 
to instabilities and partial merger of the bubbles, a complex connected structure develops. The 
parametrized DDT criterion is met first at 0.724 s after ignition at the outer edge of the flame front. 
Here, the detonation is triggered by initiating the corresponding level set, as shown in the left panel 
of Fig. ^. The fact that the outer edge of the deflagration flame is selected by the imposed DDT 
criterion is not surprising. Turbulence is generated preferentially at large buoyant bubbles and at the 
same time the density is lowest at the outermost parts of the flame making the flame broadest here. 
This favors the transition to the distributed burning regime and thus the parametrized initiation of 
the detonation wave. 

The center and right panels of Fig. ^ show the propagation of the detonation wave. Since it 
cannot cross ash regions, it wraps around the corrugated deflagration structure burning towards the 
star's center. In this way, it takes about 0.2 s before it arrives at the center of the WD. Meanwhile, 
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the star keeps expanding and the deflagration continues in regions not yet reached by the detonation 
wave. Consequently, the density of the fuel ahead of the detonation drops quickly after passing 
the WD's center. This way, burning stalls shortly before the detonation reaches the far side of 
the deflagration structure. Nonetheless, this still implies burning of most of the WD, since these 
deflagration features have already reached the low density edge of the star. 

Such simulations emphasize that the propagation of the detonation wave has to compete with 
the expansion of the WD caused by the previous deflagration phase. The outcome of this compe- 
tition with respect to the completeness of the burning, the ^^Ni production (an thus the brightness 
of the resulting event), as well as the energy release depends strongly on the efficiency of burning 
in the deflagration phase. This can be modified by varying the flame ignition geometry [53] and 
leads to models that are capable of covering the global characteristics of normal to bright observed 
SNe la [53]. Weak deflagrations lead to less expansion of the WD before the detonation is triggered 
and the detonation thus encounters a large amount of unburnt material at relatively high densities. 
Therefore it significantly contributes to the overall burning leading to very bright and energetic 
events dominated by the detonation phase. In strong deflagrations, in contrast, the high-density 
burning takes place almost exclusively in the deflagration phase. The detonation thus contributes 
little to the mass of radioactive nickel resulting in a rather dim and moderately energetic explosion 
which bears the imprint of the deflagration phase. Interestingly, however, the detonation at lower 
densities leads to a layer of intermediate mass elements surrounding the iron-group rich ejecta. 
Changing the flame ignition (for instance by assuming a multi-spot ignition scenario [25, 26]) it is 
possible to smoothly shift the characteristics of the result from a deflagration-type explosion to a 
detonation-dominated event. 

6. Conclusions 

As a conclusion, a speculation on the overall picture of thermonuclear supernovae as a model 
of SNe la can be provided based on [93]. Spectra of a set of well-observed SNe la allowed to derive 
the composition of the ejecta of the events. The result is shown in Figure ^ and indicates that the 
mass of stable iron group elements is roughly constant for all events (~0.2Mo). The mass of 
radioactive ^^Ni determines (as expected) the brightness, and the total mass of iron group elements 
shows a linear correlation with brighness. Interestingly, in all supemovae roughly the same total 
mass of burnt material is found, i.e. the dimmer events compensate less iron group elements in the 
ejecta by a larger mass of intermediate-mass elements. 

In the theoretical picture outlined above, a possible interpretation is that the dimmer exam- 
ples of the normal events (corresponding to the central part of the diagram in Figure around 
Ami5(S) ~ 1.5 mag) could be accounted for by pure deflagrations. However, because of the diffi- 
culty to synthesize large amounts of intermediate mass elements in the outer layers of the ejecta 
in deflagration models, an alternative interpretation in the context of delayed detonations seems 
more realistic. As discussed in Sect. ^^ , changing the flame ignition configuration these models 
facilitates a smooth shift from the characteristics of the model being dominated by the deflagration 
phase to a detonation-dominated event. In this picture, the SNe la in the central part of the diagram 
in Figure ^ could be explained by the weaker and dimmer models with emphasis on the deflagration 
phase while the detonation-dominated models seem promising for reproducing the bright and ener- 
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Figure 7: Distribution of the principal isotopic groups in SNe la. The enclosed mass (linked to velocity 
via the W7 explosion model) of different burning products is shown versus decline-rate parameter Ami^iB) 
(a proxy for SN luminosity). The mass of stable ^'^Fe+^^Ni for each SN is indicated in black; that of ^^Ni 
is shown in purple, and the sum of these (total NSE mass) is indicated in orange colors. Turquoise crosses 
show the sum of NSE and IME mass, indicating the total mass burned. The IME mass is given in turquoise 
color (from [93]). 



getic events (on the left-hand side of the plot in Figure 0). In this hypothetic picture, however, the 
peculiar sub-luminous SNe la located to the right in the plot remain unexplained. Further modeling 
efforts and a detailed derivation of synthetic observables from the results of numerical simulations 
will provide a way of testing this conjectured scenario. Alternative thermonuclear supernova mod- 
els, such as gravitationally confined detonations and pulsational reverse detonations, are examined 
on the basis of comparison with observations as well. While the former scenario seems to generi- 
cally lead to very bright events only, a potential problem of the latter model is the large amount of 
iron-group elements in the outer layer of the ejecta [94]. Clearly, further exploration seems neces- 
sary to settle the question whether these models can reproduce normal SNe la or lead to peculiar 
events. 
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